Quench dynamics of correlated quantum dots 
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We study the relaxation dynamics of a quantum dot with local Coulomb correlations coupled 
to two noninteracting leads which are held in grandcanonical equilibrium. Only charge degrees of 
freedom are considered and the dot is described by a model which in the scaling limit becomes 
equivalent to the interacting resonant level model. The time evolution of the current and dot 
occupancy resulting out of changes of the dot-lead coupling, the dots onsite energy, or the charging 
energy are studied. Abrupt and smooth parameter changes as well as setups with and without 
driving bias voltage are considered. For biased dots we investigate the often studied response after 
turning on the dot-lead coupling but also the experimentally more relevant case in which the voltage 
is turned on. We identify and explain a variety of interesting many-body effects and clarify the role 
of initial correlations. 

PACS numbers: 05.10.Cc, 05.60.Gg, 72.10.Fk, 73.63.Kv 



I. INTRODUCTION 

In condensed matter systems at low temperatures two- 
particle interactions lead to a variety of interesting quan- 
tum many-body effects which are not accessible by per- 
turbative methodsj^ While in equilibrium one already 
has an elaborate understanding of many aspects of corre- 
lations in the field of nonequilibrium physics this stage is 
not reached yet. One can distinguish between the physics 
of the possible (nonequilibrium) steady state and that of 
the dynamics. Here we are interested in the time evolu- 
tion resulting out of parameter changes in quantum dot 
systems with a few states tunnel coupled to two noninter- 
acting leads. The dot electrons interact via the charging 
energy and the leads are held in grandcanonical equi- 
librium; possibly with two different chemical potentials, 
that is an applied bias voltage. The parameter changes 
might be abrupt or smooth on a scale set by the system's 
inherent time scales, with the ultimate limit of adiabatic 
changes; abrupt changes are often denoted as quenches. 
The progress in nanostructuring techniques allows to ex- 
perimentally investigate quantum-dot setups with time 
dependent parameters in a rather controlled way— They 
are of interest due to their potential for conventional as 
well as quantum information processing. 

Several many-body methods successfully used in equi- 
librium were recently extended allowing to study the dy- 
namics of correlated quantum dots out of nonequilibrium 
states; some of those being (mainly) analytical*- 7 others 
purely numerical.—^— Only a few of these techniques 
have been used to treat Hamiltonians with explicitely 
time dependent parameters. We here apply a recently 
developed approximate semi-analytical method 13 which 
is based on an extension of the functional renormalization 
group (RG) approach^ to Keldysh Green functions^ - — 
It allows to study transport through quantum dots with 
a few energy levels, arbitrary time dependence of the 
single-particle parameters as well as the two-particle in- 
teraction U. The approximation of Ref. [l3| is particu- 



larly suited for dots which only have charge degrees of 
freedom, but extensions towards spinful systems showing 
the Kondo effect are envisaged. In equilibrium and for 
the bias voltage driven nonequilibrium steady state func- 
tional RG based schemes which capture aspects of the 
(spin) Kondo effect have already been developed.— ~— 
We here restrict ourselves to spinless systems and con- 
sider the interacting resonant level model (IRLM) 18 ' 23 " — 
and variants of it. 

For small to intermediate \U\ our method was shown 
to provide reliable results for the dynamics of the IRLM 
on all time scales when applied to the standard setup of 
an initially uncorrelated and decoupled dot of arbitrary 
filling and arbitrary level position e, which at time t = 
is coupled to two grandcanonical leads held at different 
chemical potentials \il = V/2 and hr = -V/2; U is turned 
on at t = as well. Here V denotes the bias voltage. 
We use units with e = h = c = 1. The underlying ini- 
tial state will be denoted as the standard initial state in 
the following and the quench protocol (turning on the 
level- lead coupling) as protocol (i). The time evolution 
of the dot occupancy and the current is characterized by 
an exponential decay towards the (7-dependent steady 
state value with two rates, which for small \U\ are related 
by a factor of two. The decay is modulated by oscilla- 
tions with frequencies given by |e ren ± V/2| and power-law 
corrections,— where e rcn denotes the (renormalized) level 
position. The functional RG^ 3 - in particular captures the 
crucial interaction induced renormalization of the decay 
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Here we study additional protocols for the parameter 
changes leading to much richer dynamics. In protocol (ii) 
we start out in the standard initial state with e = and 
let the system evolve into its steady state which has dot 
filling 1/2. We then quench the onsite energy to a large 
value. For finite bias voltage the resulting exponential 
decay of the occupancy towards its (new) steady-state 
value is repeatedly interrupted by plateau-like structures. 
They result from the interference of the oscillatory cur- 
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rent and the oscillatory occupancies of the sites next to 
the central dot one. For U = the plateaus vanish. Pro- 
tocol (iii) is related to (ii) with the crucial difference that 
the onsite energy is raised gradually over an interval T q . 
For V = and not too small T q the occupancy smoothly 
decreases and oscillations are damped; the system adia- 
batically follows the parameter change. No characteristic 
correlation effects appear. This changes for V > 0. For 
sufficiently large T q we observe an initial increase of the 
dot occupancy when e is gradually raised which vanishes 
for U -> 0. Studying the renormalization of the dot-lead 
couplings we obtain a detailed understanding of this sur- 
prising interaction effect. For both protocols we clarify 
the role of initial correlations^- by turning on U not at 
t = but at the time at which e is raised. We find that in 
the present setup correlations build up rather quickly. In 
protocol (iv) we again start in the standard initial state 
with large e as well as e = and V = 0. After the system 
reached the steady state we abruptly turn on a bias volt- 
age V > 0. To allow for this experimentally more relevant 
situation (as compared to protocol (i)) we have to slightly 
extend the formalism presented in Ref. [l3l Interestingly, 
for large e the initial dynamics of the dot occupancy re- 
sulting out of this quench differs sizably from the one of 
protocol (i) with the bias being present at t = 0. This 
already holds for {7 = 0. For e = we identify interesting 
correlation effects. The protocols studied here exemplify 
the exceptional flexibility of our approach to the quench 
dynamics of correlated quantum dots. Others can be 
easily investigated as well. 

Our paper is organized as follows. In Sect. [IT] we intro- 
duce the quantum dot model and present the extension 
of the functional RG approach to allow for varying the 
bias voltage with time. Our results for the dynamics of 
the different quench protocols are presented in Sect. IIIII 
We conclude with a brief summary in Sect. [TV] 

II. MODEL AND METHOD 
A. The quantum dot model 

Our model is given by the Hamiltonian 

H(t) = H dot (t)+ Y, [#i cad + 7C up (i)L (1) 

a=L,R 

The dot part H dot (t) consists of three lattice sites with 

H do \t) = e{t)n 2 -U{t){^ + n 2 + ^ 

+ T(t)(d\d 2 + dld 3 + R.c.), (2) 
H iat (t) = U(t)( ni n 2 +n 2 n 3 ), (3) 

in standard second quantized notation. Here rij = djdj 
denotes the occupancy operator of the spinless fermionic 
level j. The levels (sites) are connected locally through a 
hopping matrix element r > and a Coulomb interaction 



U > 0. Only the central site 2 can be moved in energy 
by changing the 'gate voltage' e. The second term in 
the single-particle part of the Hamiltonian is added for 
mere convenience so that e = corresponds to the point 
of particle-hole symmetry (half filling of the central site 
2 in equilibrium). The two leads a = L, R are modeled as 
noninteracting, 

ffi cad = E^4^ Q . (4) 

The left lead is tunnel-coupled to side 1 and the right one 
to site 3 through 

C up W=i«WEi^ + H.c., (5) 

with j£ = 1 and jn = 3. Wc explicitly allow for a time 
dependence of the parameters of H dot and H^ up . Wc 
are here interested in the parameter regime in which 
band effects do not matter and consequently assume 
the leads to be in the wide band limit (for details see 
Ref. liH). The level- lead coupling is then parameterized 
by T a (t',t) = nt^(t')t a (t)D a , with the constant density 
of states D a of lead a. To reduce the number of param- 
eters we focus on = Tr = T. The leads are held at 
grandcanonical equilibrium with temperature T = and 
chemical potential Hl/r = ±V/2, with the bias voltage 
v X). 

In the scaling limit t, |e], \U\, V « T the model is equiv- 
alent to the field theoretical IRLM^ 3 - with a single energy 
level as the first and third dot sites can effectively be in- 
corporated into the leads . 13 ' 18 Here we do not solely focus 
on the scaling limit but also consider setups with t % T 
and U £ T. For r £ T our model describes a dot with 
three levels. The linear geometry and the coupling of 
sites 1 and 3 to the leads implies a specific coupling of 
these three energy levels. 

In Ref. [Ill we gave a detailed account on how to ob- 
tain approximate expressions for the retarded, advanced 
and Kcldysh Green functions and the corresponding self- 
energies for Hamiltonians of the above type using func- 
tional RG. From the Green functions and self-energies 
the occupancies hj of the dot sites and the left and right 
current J a can be computed using standard formulas also 
given in Ref. [T3 . Consequently we here refrain from pre- 
senting any technical details on the functional RG and 
refer the interested reader to Ref. Il3l . Our approximation 
consists in keeping only the lowest order diagrams in the 
two-particle interaction U on the right hand side of the 
RG flow equations for the vertex functions. Solving these 
differential equations however implies resummations of 
classes of diagrams leading to results beyond simple per- 
turbation theory. We e.g. capture the renormalization 
of the decay rates characterized by a power-law behav- 
ior in the bare rates with a {7-dependent exponen t 18 ' 23 
as well as the power-law correction to the exponential 
decay in time of n 2 (t) within protocol (i) with an in- 
teraction dependent exponent . 13 ' 26 Both exponents come 
out correctly to leading order within our approximation 
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and the method is thus reliable for \U\ £ T. We em- 
phasize, that within the functional RG based scheme we 
can study the dynamics on arbitrary time scales. This 
has to be contrasted to numerical approaches to time- 
dependent nonequilibrium transport through correlated 
quantum dots with which large times can currently not 
be reached. A detailed comparison of results for the stan- 
dard protocol (i) to those obtained by other approaches 
can be found in Ref. HH For completeness and compari- 
son the physics of this protocol is briefly summarized in 
Sect. IIII A"l Within our approach the two-particle interac- 
tion leads to a (time dependent) renormalization of the 
single-particle parameters. To compute observables at 
the end of the RG flow one has to solve a single-particle 
problem with renormalized parameters indicated by the 
superscript 'ren'. 

In protocol (iv) the bias voltage is turned on at a time 
t q > 0. This type of time dependence is not directly 
captured by the formalism of Ref. [TH which we thus have 
to supplement. The necessary steps are described in the 
next subsection. 



B. Turning on the voltage 

In Ref. [H we discussed how to handle time depen- 
dent couplings to the reservoirs via the time dependent 
hybridization T a (t',t). We here show that a gauge ar- 
gument can be used to shift time dependent chemical 
potentials /i Q (£) into time dependent hybridizations; this 
way protocol (iv) is accessible to the approach of Ref. [l3| 
as well. 

For structureless wide band limit leads considered here 
raising the lead filling, and thus the chemical potential, 
is completely equivalent to raising all the energy levels of 
the leads. Therefore, we can replace the lead part of the 
Hamiltonian Eq. Q by 



H, 



lead 



•Ma(*)) 



(6) 



and take the leads distribution functions as Fermi func- 
tions with zero chemical potential. We then apply the 
gauge transformation 



G(t) 



(7) 



where fl a (t) is a suitable gauge field with ^(O) = 0. The 
transformed Hamiltonian reads 



H(t) = G{t)H{t)G\t) + iG{t)G\t), 



(8) 



which shows that for setting Cl a (t) = ^, a {t) the addi- 
tional term in Eq. (|6]) originating from the time de- 
pendent chemical potential cancels with iG(t)G'(t) = 
-T a T ko n a (t)c{ a C ka . Since [H™ up (t), G(t)] * the 
hopping amplitude describing the level-lead coupling in 
H acquires an additional phase factor 
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FIG. 1: (Color online) Time evolution of the occupancy 
and the current leaving the left reservoir for the IRLM in 
protocol (i). The parameters are t/T = 0.0025, e/Tk = 10, and 
V/Tk = 10. The universal energy scale Tk is renormalized by 
the local Coulomb interaction U (see the text). 



where f2 a (0 = Jq dt'fi a (t'). The transformed Hamilto- 
nian can then be treated as outlined in Ref. Il3l 



III. RESULTS FOR THE QUENCH DYNAMICS 

In this section we present and discuss our functional 
RG results for the time evolution of the occupancies nj 
and the currents J a induced by the quench protocols de- 
fined in the introduction. The focus is on the correla- 
tion effects. For comparison and to gain a deeper un- 
derstanding we present U = results when necessary. 
Already now we emphasize that some of the interaction 
effects cannot directly be linked to the renormalization 
of the decay rate, which in the standard setup proto- 
col (i) was found to be the most dominant effect of the 
two-particle interaction . 18 i 24 ~ — To keep our presentation 
self-contained the physics of protocol (i) is briefly sum- 
marized in the following. 



A. Protocol (i) 

In protocol (i) we start out with an initially empty 
dot decoupled from the two grandcanonical leads held at 
chemical potentials ^l/r = ±V/2. The coupling to the 
reservoirs is turned on at time 0. Figure [T] shows the 
time evolution of the level occupancy fii and the current 
Jl leaving the left lead in the scaling limit (the IRLM) 
for t/T = 0.0025, e/T K = 10, and V/T K = 10. Here T K 
is the universal energy scale defined by the equilibrium 
charge susceptibility-i 8 -^^^ 



dn 2 
~~d7 



ttTk 



(10) 



t a {t) ^ t a {ty^ l \ 



(9) 



For U = it is given by 4r 2 /r, but gets strongly renor- 
malized by the interaction . 18 ' 23 " — It increases for small 
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to intermediate repulsive interactions. The exponential 
decay of n 2 and J a towards their respective [/-dependent 
steady-state values is characterized by two relaxation 
rates. For small U they are approximately given by Tk 
and Tk/2. The decay with the smaller rate is accompa- 
nied by a power-law correction in t with a [/-dependent 



exponent 



13,18.26 



Since this is only a small correction to 



the exponential decay the main effect of the interaction is 
the strong enhancement of the relaxation rates for U > 0. 
For |e ron ± V/2\ » T K ,l/t the decay of the occupancy 
7T, 2 (i) is overlayed by oscillations with the two frequen- 
cies |e rcn ± V/2\, where the renormalization of the on- 
site energy e -*■ e rcn is of order (U/T) 2 and thus small. 
In the current leaving the left (right) lead the frequency 
|e ren -F/2| {\e rcn + V/2\) dominates Comparing the 

curves for U = and U > in Fig. [TJ it is apparent that 
the two-particle interaction leads to an increase of the 
oscillation amplitudes. 



B. Protocol (ii) 

We next study the response of a system, which after 
being prepared according to protocol (i) has reached a 
steady state, to the sudden increase of e from e; n ;/r = 
to the rather large value efi na i/r = 2 at time t q > 0. We 
choose r/T = 0.5, which is no longer deep in the scaling 
limit. The time-evolution for t < t q (see Fig. [2]) still shows 
the main features found for the IRLM discussed in the 
last subsection. 

Figure [5] shows the time evolution of the occupancies of 
the three sites as well as Jj, for V/T = 2. The occupancy 
of the middle site and the current drop quickly after the e- 
quench at Tt q = 20 as the onsite energy is raised far above 
the chemical potential. Whereas in the noninteracting 
case fi2(i) smoothly decreases (up to weakly developed 
shoulders; dashed line in the upper panel) this depletion 
is periodically interrupted for U > (solid line; see the 
ellipses in the upper panel). This effect, which clearly 
goes beyond the renormalization of the decay rates, can 
be understood as follows. The interaction enhances the 
coherence of the system leading to an increase of the oscil- 
lation amplitudes as discussed for protocol (i). Plateau- 
like features show up for time regimes where maxima of 
the current leaving the left reservoir coincide with the de- 
creasing part of the oscillatory occupancy of the first site 
n\(t). The combination of these tendencies temporarily 
stops the depletion process on the central site leading to 
the plateaus. These time regimes are marked with A in 
Fig. [3J For times in which minima in the left lead's cur- 
rent coincide with an increase of fix, marked with B in 
Fig. [5J the depiction of the central site is even enhanced 
compared to the noninteracting case. The comparison to 
the curve for U/T = 1 but V = (dashed-dotted line in 
the upper panel) shows that a finite bias voltage is vital 
for the development of the plateaus. 

To understand the periodic interception of the decay 
on a more quantitative level we determine the frequencies 
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FIG. 2: (Color online) Top: Time evolution of the occupancy 
n-2 in protocol (ii). The quench from ei n i/T = to £fl na i/r = 2 
is performed at time Tt q = 20. At this time the system has 
reached a steady state after initially being prepared accord- 
ing to protocol (i). The other parameters are t/T = 0.5 and 
V/T = 2. The plateau-like features interrupting the exponen- 
tial decay beyond Tt q = 20 are highlighted by ellipses. For 
comparison n,2{t) for U/T = 1 and V = is shown as well. 
Bottom: Time evolution of n 1 / 3 and Jl- 



expected for t > t q from protocol (i). For the noninter- 
acting case those are given by Eq. (A. 5) of Ref. [HI con- 
sidering the poles associated to the smallest relaxation 
rates (which dominate at large times). As a first ap- 
proximation one can use the same formula with the bare 
parameters replaced by the renormalized ones taken for 
t -> oo. This gives an estimate of the corresponding fre- 
quencies in the interacting case. For the parameters of 
Fig. [J] we obtain 



w rcn + V/2 « 3.23747r 
u ien -V/2* 1.23747T. 



(11) 
(12) 



These two frequencies correspond roughly to the frequen- 
cies associated to the periodic appearance of the plateaus 
(smaller frequency) as well as the oscillation frequencies 
within the plateaus (higher frequency). 



C. Protocol (iii) 

We now focus on the case in which the change of 
the onsite energy described in the last subsection is not 
abrupt, but is performed with constant velocity over the 



time interval T a 



centered around t q from 



i/r = to 



efiniai/r = 2. To ensure that the system has fully reached 
its steady state at time t = t q - T q /2 we modify the initial 
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FIG. 3: (Color online) Time evolution of n2 for a change of 
e performed with constant velocity in the time interval T q 
(protocol (hi)). The parameter are chosen as in Fig. [2] but 
with V = 0. Arrows above the graph indicate the time interval 
over which the change of e is performed. 
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FIG. 4: (Color online) The same as in Fig.[3]but with V/T = 2. 
For comparison the time evolution after an abrupt quench (see 
Fig. [2} is shown as well. The inset shows a zoom-in of the 
region in which n2 (initially) increases. Arrows above the 
graph indicate the time interval over which the change of e is 
performed. 



preparation according to protocol (i) by starting at t = 
with an initial density matrix with occupancy 1/2 on its 
diagonal. 

We start the discussion with the unbiased case V = 
and choose the other parameters as r/T = 0.5 and U/T = 
1. Figure [3] shows the time evolution of t%2 for different T q 
(see the arrows above the graph) with Tt q = 30 fixed. As 



1.5 



T 



OV;n =12.12 

12 q 

l<"rV;rT =12.12 

23 q 




10 



20 



30 

rt 



40 



50 



60 



FIG. 5: (Color online) Renormalized hopping amplitudes 
across bonds (1,2) and (2,3) as a function of time (normal- 
ized to the bare time-independet value). The parameters are 
chosen as in Fig. |4]with TT q = 12.12. 



in Fig.|4j Interestingly, for U > and sufficiently large T q , 
77,2 (i) increases for times slightly larger than t q -T q /2 (see 
the dashed and dashed-dotted lines). For larger T q the 
effect becomes more pronounced. There are two compet- 
ing processes determining the occupancy n 2 (t). (a) The 
increase of the onsite energy favors a depletion process 
and dominates in the case where the quench is abrupt, 
(b) At finite bias voltage the interaction induced renor- 
malization of the hopping amplitudes across bonds (1,2) 
and (2,3) is not equal as shown in Fig. [5j The hopping 
amplitude across bond (1,2) (solid line), which is closer to 
the lead held at higher chemical potential, increases more 
than the hopping across (2,3) (dashed line), by which 
charge is transported to the lead held at lower chemical 
potential. This induces a tendency towards an increase of 
the charge on the central site. When the increase of e is 
performed slow enough the increase of 77,2 due to the hop- 
ping asymmetry can initially overcompensate the deple- 
tion. At larger times process (a) wins and the occupancy 
decreases exponentially (with overlayed oscillations). We 
emphasize that the (initial) increase results from the in- 
terplay of the finite bias voltage V > 0, the two-particle 
interaction U > 0, and the parameter change not being 
abrupt. It vanishes when either of these conditions is not 
fulfilled. 



D. The role of initial correlations 

In protocols (ii) and (iii) the two-particle interaction is 
turned on at t = 0. Correlations are thus fully developed 
at the time of the abrupt or gradual increase of e. To 



expected the amplitude of the oscillations decreases with investigate the role of these initial (with respect to t q or 
increasing T q . In the limit T q -+ 00 the system adiabat- T q - t q /2) correlations in the quench dynamics we now 



ically follows the parameter change and the oscillations 
disappear. 

Next we consider the biased case with V/T = 2 shown 



modify the 'initially correlated' protocols (ii) and (iii) 
and abruptly turn on the interaction at t q or T q - t q /2, 
respectively. We will refer to these modified protocols as 
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FIG. 6: (Color online) Absolute value of the difference of the 
derivatives of fi2{t) for the initially correlated and uncorre- 
cted protocol (ii). The parameters are as in Fig. [2] The 
dashed line displays an exponential decay with the rate set 
by the final Hamiltonian. Note the linear-logarithmic scale. 
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FIG. 7: (Color online) Comparison between the time evolu- 
tion of U2 for protocol (iii) with and without initial correla- 
tions. A zoom-in of the interesting time regime in which e is 
raised is shown. 

'initially uncorrelated'. 

We start out with the comparison of the time evolution 
of the initially correlated and uncorrelated protocol (ii) . 
It is not surprising that for our tunnel coupled quantum 
dot setup the occupancies n 2 (t) of the two situations be- 
come equal exponentially fast. In particular, the plateau- 
like features discussed in connection with Fig. [5] are also 
present in the initially uncorrelated setup. To make the 
differences clearly visible in Fig. [5] we show the absolute 
value of the derivative of the difference of the two occu- 
pancies fii for t > t q on a linear-logarithmic scale. The 
envelope shows a linear behavior which proofs that the 
correlations build up on an exponential scale. In fact, 
this scale is found to be the decay rate set by the final 
Hamiltonian (after the quench). An exponential decay 
with this rate is shown as the dashed line in Fig. [BJ 

In protocol (iii) we found that the correlation induced 
increase of sets in right after the time t q - T q /2 at 

which e starts to be gradually raised. It is thus rea- 



sonable to expect more obvious differences in the dy- 
namics of the central sites occupancy for the initially 
correlated and uncorrelated states of this protocol. As 
Fig. [7] shows the local maximum described in subsection 
IIII CI decreases in height if initial correlations are absent 
(dashed-dotted line with squares) up to the extend that 
the maximum disappears for the faster parameter change 
with TT q = 2.52 (dashed line with squares). Additionally, 
the short time behavior right after the parameter change 
has been initiated changes significantly. In the slowest 
increase of e with FT q = 12.12 (dashed-dotted line with 
squares) first a minimum develops. The initial decrease 
of fii reflects the crucial importance of two-particle cor- 
relations for the accumulation of charge on the central 
dot. After significant correlations have build up (in time) 
n<i(i) increases as found for the initially correlated proto- 
col (iii) . The long-time behavior of the initially correlated 
and uncorrelated protocols (iii) is indistinguishable. 



E. Protocol (iv) 

When studying time-dependent currents and level oc- 
cupancies of correlated quantum dots in most theoretical 
studies the response towards turning on the level-lead 
couplings is considered. In experiments on the other 
hand one generically starts out of an equilibrium state 
in which the leads are coupled and at a certain time 
turns on the driving bias voltage. Next we therefore con- 
sider this experimentally motivated protocol (iv) in the 
scaling limit (the IRLM). For this we start out in the 
standard initial state with V = and empty dot follow- 
ing protocol (i). After a steady state is reached the bias 
voltage is quenched to a positive value. In Fig. |5]n2(t) 
for the abrupt voltage quench is compared to the evolu- 
tion out of the standard initial state (with V > 0) as dis- 
cussed in subsection llll Al The parameters are V = 10Tk, 
r/r = 0.0025, e = 10T K , and U/T = 0.2 as well as 0. The 
occupancy for the two cases shows sizable differences. 
While in the evolution out of the standard initial state 
the two frequencies \e Ten ±V/2\ dominate (see subsec- 
tion IIII A[) the shape is deformed towards a (damped) 
triangular-like time dependence when the bias is only 
turned on at a later time; instead of two, many frequen- 
cies uj with e rcn - V/2 < \w\ < e ren + V/2 contribute to 
this line shape. The deformation is already apparent for 
U = 0. Furthermore, the amplitude of the oscillations 
is reduced when the bias voltage is quenched. In con- 
trast to n%(£), the line shapes of the current Jl(£) in 
protocols (i) and (iv) shown in Fig. [5] are rather similar. 
The only prominent difference is the amplitude which de- 
creases if the bias is turned on at a later time (also at 
U = 0). The correlation effects of protocol (iv) for large e 
are the renormalization of the decay rate 7k (note that 
the time axis of Figs. [8] and [9] is given in units of 7k) and 
the slightly increased amplitude of the oscillation, both 
known from protocol (i). 

In addition to the case of a large onsite energy we study 




FIG. 8: (Color online) Time evolution of the occupancy «2 
for a quench in the bias voltage (protocol (iv)) and the level- 
lead coupling (protocol (i)). The parameters are V = lOTk, 
r/r = 0.0025, e = 10T K , and U/F = 0.2 as well as 0. For a 
better comparison the data or protocol (i) are shifted along 
the time axis to the time of the voltage quench Tkt q ~ 5.5. 
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FIG. 9: (Color online) The same as in Fig. [8] but for the 
current Jl(*). 



the dynamics for e = and initial occupancy n^it = 0) = 0. 
The steady-state occupancy equals 1/2 and it is useful to 
consider the quantity \2ri2(t) - 1| on a logarithmic scale, 
which enhances the visibility of many features. For V = 
the IRLM is equivalent to the ohmic spin-boson model 28 
and |2n 2 (£) - 1| corresponds to the absolute value of the 
z-component of the spin expectation value. Figure 1101 
shows the time evolution of |2ri2(t)-l| out of the standard 
initial state and a subsequent quench in the bias voltage 
from V = to V = 10T K at T K t q = 10. The noninteracting 
case as well as U/T = 0.1 are depicted. The hopping is 
r/r = 0.0025. For U = 0, |2n 2 (i) - 1| goes to zero purely 
exponentially with a single dominant rate Tk(U = 0); 
straight line (dashed) on the linear-logarithmic scale of 
Fig. [TO] This can be understood from the explicit an- 
alytical expression for ^(i) in the noninteracting case 
e.g. given in Eq. (79) of Ref. [H. For e = the two os- 
cillatory terms cancel. For U > one observes the expo- 
nentially damped oscillations with frequency ~ U of the 
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FIG. 10: (Color online) Time evolution of \2n 2 (t) - 1| for 
a quench in the bias voltage from V = to V = 10Tk at 
Tktq = 10. Results for U/F = 0.1 as well as are shown. 
The other parameters are r/r = 0.0025 and e = 0.0. For 
comparison a graph for a setup in which the interaction is 
turned on simultaneously with the bias is shown as well. Note 
the linear-logarithmic scale. The time axis is given in units of 
T K (U = 0) for the dashed line and T K (U = OAF) for the solid 
and the dashed-dotted line. This leads to the two different 
slopes and the kink close to t q in the dashed-dotted line. 



coherent phase of the spin-boson mode l ' ' before the 
bias quench (dip for t < < q of the solid line) and additional 
oscillations with much higher frequencies after this. Ad- 
ditionally, the dynamics of a setup in which V as well as 
U are turned on at t q (no correlations for t < f q ) is shown 
as the dashed-dotted line. Although the relaxation rate 
appears to change to its renormalized value almost in- 
stantly as the interaction is turned on (kink at t « t q and 
the same slope as the other lines for t>t q ; note that the 
time is given in units of the interacting 7k for all times), 
the dynamics differs significantly compared to the one 
of the solid line; on the scale of the plot no oscillations 
are visible. We conclude that in this example the build- 
up of the correlations quickly changes the relaxation rate 
but does not result in significant oscillations as found for 
the initially correlated case. We here concentrated on 
the occupancy, as the dynamics of Jl proves less inter- 
esting and is characterized by an exponential approach 
of the steady-state value modulated by oscillations with 
frequency V/2. 



IV. SUMMARY 

We have studied the dynamics of the dot occupancy 
and the current resulting out of abrupt quenches and 
smooth parameter changes in the IRLM and variants of 
this. The time evolution shows characteristics which de- 
pend on the protocol with which the different parameters 
are varied in time. In particular, correlations can lead to 
a rich dynamics showing interesting effects which go be- 
yond the renormalization of the decay rate found in the 
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standard protocol (i). Using our approximate functional 
RG approach we were able to identify and explain sev- 
eral of those. Our method is very flexible, can be used for 
all time scales, and other time dependences of the model 
parameters than the ones considered here can be treated 
as well. We exemplified the role of initial correlations 
and determined the difference between the dynamics re- 
sulting out of switching-on the level-lead coupling in the 
presence of a finite bias voltage and the experimentally 
more relevant case of turning-on the bias voltage in an 



equilibrium state with coupled leads. 
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